Measurement of the neutron F 2 structure function via spectator tagging with CLAS 
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We report on the first measurement of the F2 structure function of the neutron from the semi- 
inclusive scattering of electrons from deuterium, with low-momentum protons detected in the back- 
ward hemisphere. Restricting the momentum of the spectator protons to < 100 MeV/c and their 
angles to > 100° relative to the momentum transfer allows an interpretation of the process in terms 
of scattering from nearly on-shell neutrons. The F£ data collected cover the nucleon-resonance and 
deep-inelastic regions over a wide range of Bjorken x for 0.65 < Q 2 < 4.52 GeV 2 , with uncertainties 
from nuclear corrections estimated to be less than a few percent. These measurements provide the 
first determination of the neutron to proton structure function ratio F£ /F$ at 0.2 < x < 0.8 with 
little uncertainty due to nuclear effects. 

PACS numbers: 13.60-r, 13.60.Hb, 14.20.Dh 



Structure functions of the nucleon reflect the defining 
features of QCD: asymptotic freedom at short distances 
and quark confinement at long distance scales. After four 
decades of deep-inelastic lepton scattering (DIS) mea- 
surements at facilities worldwide, an impressive quan- 
tity of data has been collected, extending over several 
orders of magnitude in Bjorken x (the fraction of the nu- 
cleon's momentum carried by the struck quark) and Q 2 
(the squared 4- momentum transfer). These data have 
provided strong constraints on the quark and gluon (or 
parton) momentum distribution functions (PDFs) of the 
nucleon. 

Although the structure of the proton has been well de- 
termined, the absence of high density, free neutron tar- 
gets has meant that neutron structure functions must 
be inferred from experiments on nuclear targets, partic- 
ularly deuterium. In regions of kinematics where most of 
the momentum resides with a single quark, x > 0.5, un- 
certainties in the nuclear corrections in deuterium result 
in large uncertainties in the extracted neutron structure 
functions [IHE]- 

Determining structure functions and PDFs at large x 
is important for several reasons. For example, one of the 
long-standing puzzles in hadronic physics is the behavior 
of the ratio of d to u quark PDFs in the proton in the 
limit x — > 1 [2|. A number of predictions have been 
made for the d/u ratio in this limit, from perturbativc 
and nonperturbative QCD arguments Q , but because of 
the lack of neutron data these have never been verified. 

A better knowledge of neutron structure functions in 
the resonance region is also needed to help unravel the full 
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isospin structure of the resonant and nonresonant contri- 
butions to the cross section, as well as to provide critical 
input for interpreting inclusive polarization asymmetry 
measurements. An important question in the resonance 
region is whether Bloom-Gilman duality holds as well for 
the neutron as it does for the proton [8, 9]. Furthermore, 
model-independent determinations of F2 arc essential for 
evaluating the efficacy of model-dependent extractions 
[Toj of F£ in the resonance region from inclusive deu- 
terium data. 

It has been suggested [ilTflli ] that one can greatly 
reduce the nuclear model uncertainties associated with 
scattering on the dcuteron by selecting events with low- 
momentum protons produced at backward kinematics 
relative to the momentum transfer. Tagging backward- 
moving spectator protons minimizes final-state interac- 
tions [141 1 1 51 ] , and the restriction to low momenta ensures 
that the scattering takes place on a nearly on-shell neu- 
tron. Furthermore, by measuring the momentum of the 
recoiling proton, one can correct for the initial motion of 
the struck neutron, all but eliminating Fermi smearing 
effects. 

In this Letter we report on the first direct extraction 
of the neutron F£ structure function by tagging specta- 
tor protons in semi-inclusive electron scattering from the 
dcuteron. In the impulse approximation, where the vir- 
tual photon scatters incoherently from a single nucleon, 
the differential cross section for the reaction d(e,e'p s )X 
is given by [H, EH 

da _ 2a 2 f x 2 y 2 M 2 \ 

dxdQ 2 d 3 Ps /E s xQ 4 V ~ V ~ Q 2 J 

x(V 2 d + 2tan 2 ^-^), (1) 
where a is the fine structure constant, p s = \p s \ and 
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E s = \J M 2 + p 2 are the spectator nuclcon momentum 
and energy in the laboratory frame, and M is the nu- 
cleon mass. Here x — Q 2 /2Mv is the Bjorken scaling 
variable, with v the energy transfer to the deuteron, and 



Q 



-q is the square of the exchanged virtual photon 



4- momentum vector q. The variable y = v/E is the frac- 
tion of the incident electron energy E transferred, and 
9 is the electron scattering angle. Additional structure 
functions that vanish after integration over the azimuthal 
angle of the spectator have been omitted in Eq. ([TJ . 

The semi- inclusive deuteron structure functions Tf 2 
are in general functions of four variables, J~f 2 = 
J-f 2 (x, Q 2 1 a s,P't), where a s = (E s —p z s )/M is the frac- 
tion of the deuteron's light-cone momentum carried by 
the spectator proton, and p z and pjr are its longitudinal 
and transverse momenta, respectively. In the impulse ap- 
proximation the functions Tf 2 are related to the (effec- 
tive) neutron structure functions F™ 2 and the deuteron 
spectral function S(a.s, jpj~); in the limit of large Q 2 and 
small pjr /M one has [l2| 



F«<xS(a s ,pi)Fr 2 (x*,Q 2 ,p 2 ), 



(2) 



where x* = Q 2 /2p ■ q w x/(2 — a s ) is the Bjorken scaling 
variable of the struck neutron in the deuteron, and p 2 = 
(Aid — Eg) 2 — p 2 is its virtuality, with the deuteron 
mass. The spectral function is proportional to the square 
of the deuteron wave function. In terms of a;* the inferred 
invariant mass squared of the struck neutron remnant is 
given by W* 2 = (p+q) 2 = p 2 + Q 2 (l-x*)/x* , in contrast 
to the usual definition of W 2 — M 2 + Q 2 (l —x)/x for a 
free nucleon. 

For inclusive scattering on the deuteron one integrates 
Eq. {TJ) over all spectator momenta p s and expresses the 
extracted structure function in terms of the variables x 
or W; for the tagged reaction the detection of a proton 
at specific kinematics selects a fixed x* and W* . More- 
over, the restriction to backward-moving protons serves 
to minimize the probability of the recoil proton rescat- 
tcring with the debris of the struck neutron. Calcula- 
tions within hadronization models suggest [3 EH that 
for spectator momenta below ~ 100 MeV/c final-state 
interaction effects distort the spectral function by < 5%, 
provided that spectator angles pq are above 100°. Back- 
ward kinematics also suppresses contributions from low- 
momentum protons emanating from the hadronic debris 
of the struck neutron, which distort the spectral function 
at the < 1% level [l3|. These theoretical calculations are 
corroborated by both existing data [l6| and by our own 
analysis of the full data set [17 1 . 

Because the neutron is bound inside the deuterium nu- 
cleus with binding energy Ed = —2.2 MeV, it can never 
be exactly on-shell since p 2 — M 2 sa 2Msd — 2p 2 < 0, 
even when it is at rest. The dependence on the neutron's 
virtuality may introduce additional differences between 
the effective neutron structure functions in Eq. @ and 
their on-shell values. However, since the bound neutron 
is « 13 MeV away from its mass-shell for p s = 100 MeV/c 



(and only 7.5 MeV for p s = 70 MeV/c) the uncertainty 
introduced in extrapolating to the on-shell point is min- 
imal. Indeed, quantitative estimates of the off-shell de- 
pendence of the neutron structure functions in relativistic 
quark-spectator diquark models [l8l . [l9| and models that 
consider the effects of evaluating the structure function 
at a shifted energy transfer [2(| give corrections to the 
on-shell structure functions of < 1% forp s < 100 MeV/c. 

The BoNuS (Barely off-shell Nucleon Structure) exper- 
iment ran in 2005 using the CEBAF Large Acceptance 
Spectrometer (CLAS) [2l[ in Hall B at Jefferson Lab. 
Electrons scattered from a thin deuterium gas target were 
detected by CLAS and the spectator protons were mea- 
sured with the BoNuS Radial Time Projection Chamber 
(RTPC) [23]. Production data were taken at three beam 
energies, 2.140, 4.223 and 5.262 GeV, with an additional 
set of calibration data taken at 1.099 GeV. The kine- 
matic coverage includes final-state invariant masses from 
the quasielastic peak up to W* « 3 GeV, and momentum 
transfers Q 2 from 0.2 to « 5.0 GeV 2 . 

The RTPC reconstructed the three-dimensional tracks 
of spectator protons in a 3 cm wide annular ionization 
volume, using gaseous electron multipliers to amplify the 
ionization electrons. The signals were read out via a 
grid of conducting pads on a cylindrical outer surface 
in 114 ns increments of time, yielding up to 60 points 
in radius, azimuth and z (the distance along the beam 
direction) for each track. The 170 mm long target inside 
the 200 mm long RTPC allowed the detection of spec- 
tator protons with polar angles 20° < 6 S < 160° in the 
lab frame, covering 295° in azimuth. This provides good 
spectator acceptance over the range —0.9 < cos pq < 0.9. 
The detector was immersed in a 4 T solenoidal magnetic 
field which suppressed the electromagnetic background 
(Mollcr electrons) and bent the proton tracks. Measur- 
ing the curvature allowed the reconstruction of the proton 
momentum, and measuring the total ionization charge as- 
sociated with a track enabled the separation of protons 
from other hadrons through their specific energy loss. By 
requiring tracks to be in time with the detected electron 
(within 2 fj,s) and to trace back to the electron vertex 
in z (within 30 mm), accidental backgrounds could be 
suppressed to about 20%. Using events with a larger dis- 
tance in z between the electron and proton vertices as 
a sample of accidentals, this background was subtracted 
from the data. Details of the RTPC construction and 
performance are found in Ref. [221 ], 

The data were also corrected for pions misidentificd 
as electrons in CLAS and for electrons coming from 
pair-symmetric decays of mesons and photons. Cuts on 
y < 0.8 eliminated events with large radiative correc- 
tions. Lower limits were placed on x for each bin in 
Q 2 to remove acceptance edge effects. The low density 
of material in the path of the outgoing protons allowed 
them to be identified with momenta down to 70 MeV/c. 
The analysis was restricted to protons with momenta less 
than 100 MeV/c, and angles relative to the momentum 
transfer vector q of more than 100° - in the following 
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referred to as the kinematic bin A( VIp )p s for "very im- 
portant protons" (VIPs). 
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FIG. 1: Yield for the semi-inclusive d(e,e'p a )X reaction 
with a backward-moving spectator proton as a function of the 
invariant mass W* of the neutron debris, compared with the 
yield for the inclusive d(e, e')X reaction as a function of the 
customary kinematic variable W . Yields integrated over W 
and W* are normalized to be the same. The data are for the 
4.223 GeV beam energy and are averaged over the acceptance 
of CLAS. For backward-moving spectators W* < W , which 
explains the leftward shift of the high W* cutoff in the semi- 
inclusive spectrum with respect to the inclusive case. 

The utility of the spectator tagging method is illus- 
trated in Fig. [TJ where a typical semi-inclusive yield for 
the d(e, e'p s )X reaction is shown as a function of the in- 
variant mass W* of the neutron's hadronic debris, and 
the corresponding inclusive yield for the d(e, e')X reac- 
tion is shown as a function of the usual invariant mass 
W for a neutron struck at rest in the lab frame. The 
quasielastic and A(1232) resonance peaks are largely 
smeared out by the nuclear Fermi motion in the inclu- 
sive spectrum, whereas the neutron elastic and resonance 
peaks clearly stand out in the semi-inclusive spectrum. 
The elastic neutron peak for d(e,e'p s )X has a Gaussian 
width of 31 MeV, which is only 20% larger than that for 
a proton target measured with CLAS. 

For our final results, we formed the ratio F OX p of the 
acceptance-corrected yields for d(e, e'p s )X in the individ- 
ual W* (or x*) and Q 2 bins for a spectator proton within 
the bin A^ VIP ^p s , divided by the similarly corrected yield 
measured for d(e, e')X at the corresponding W or x, 



Re 



taggo 



d (AQ 2 , AW*, A^ Ps )/A e (Q 2 , W*) 



N incl (AQ 2 ,AW)/A e (Q 2 ,W) 



(3) 

In this ratio, the total luminosity of the experiment can- 
cels, and the corrections due to the CLAS acceptance for 
the scattered electrons, A e , largely cancel, as this enters 



the numerator and denominator at rather similar kine- 
matics. The acceptance A e was determined from the 
ratio of inclusive electron count rates and the known ed 
cross section (24|. Although A e varied by a factor of 2, 
the corrections to the ratio were less than 10% with a 3% 
uncertainty. Radiative corrections were applied to both 
the numerator and denominator based on the prescrip- 
tion by Mo and Tsai [Hj], using models 0] of F 2 ™, F 2 d 
and the ratio of longitudinal to transverse cross sections 
as input for the calculations. These also canceled to a 
large extent in the ratio and were less than 10% with a 
2% uncertainty. 

In the spectator approximation of Eq. @, the ratio 
i?exp is directly proportional to the ratio of (free) struc- 
ture functions F£ /F% multiplied by the spectral function 
S(a s ,pjr) integrated over the proton acceptance A p of 
the RTPC within the VIP cuts, 



R 



cxp 



F£(W*,Q 2 ) 
F 2 d (W,Q 2 ) 



VIP 



da s dpj A p (a s ,p^) S(a s ,p^-). 

(4) 

The integral /yip m Eq. ^ is largely independent of 



kinematics, and (F?/F£) exp = R cxp {F$ / F%) / I V1P , in 
which inland F 2 P are well-measured values parametrized 
in Rcf . [241 ■ The normalization constant iyip wa s chosen 
for the whole data set using F^/Ff = 0.695 at x = 0.3, 
where nuclear effects are small, with an uncertainty of 3% 
from the CTEQ- Jefferson Lab (CJ) global PDF fits @]. 
The rms variation in the normalization constant iyip for 
subsets in W* and Q 2 was 3.4%, which was included in 
the systematic error. The structure function {F£) exp was 
obtained by multiplying (F-f/Fj ) CX p by the values of F 2 
parametrized in Ref . [24j . The final systematic errors in- 
clude uncertainties on F 2 and F 2 and possible deviations 
from the (implicit) assumption that the longitudinal to 
transverse cross section ratios are the same for d, p and n, 
as well as residual background, acceptance and radiative 
correction uncertainties. A conservative systematic error 
of 3% was assigned to possible violations of the spectator 
assumptions due to final-state interactions and off-shell 
effects [l2T - fl5T ] . An additional 3% (rms) uncertainty arises 
from the global fit for F 2 . 

A representative sample of the neutron F 2 spectra 
is shown in Fig. [5J compared with a phenomenological 
parametrization of F 2 [24[ obtained from inclusive F d 
and F 2 data using a model of nuclear effects, and an 
extraction [l(j of F£ from recent F% and Ff data us- 



ing the nuclear smearing corrections of Ref. [25J]. (The 
complete spectra for all kinematics are published in the 
CLAS database [H[.) 

The comparison shows reasonable overall agreement 
between the BoNuS data and the model-dependent F£ 
extractions [l(| [H| from inclusive data, but highlights 
some residual discrepancies. In particular, at the lowest 
Q 2 values both the parametrization [24| and the model- 
dependent extraction [l(| underestimate the F£ data, es- 
pecially in the vicinity of the A(1232) peak. At larger Q 2 
the models are in better agreement with the data in the 
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FIG. 2: (color online). Typical F% spectra from the BoNuS 
experiment (filled circles) as a function of W* for the vari- 
ous Q 2 ranges indicated. The beam energy was 5.262 GeV 
except for the upper left plot at 4.223 GeV. For comparison 
the model-dependent extraction from inclusive F$ data (o pen 
circles) [l(| and the phenomenological model from Ref. [241 ] 
(solid curve) are also shown. The error bars on the data points 
are statistical, and the band along the abscissa represents the 
systematic error without the overall 3% normalization uncer- 
tainty or the 3% spectator approximation uncertainty. 



A region, but overestimate it somewhat in the third res- 
onance region at Q 2 ~ 2.5 GeV 2 . This suggests that 
cither the nonresonant neutron contribution assumed in 



the model [24J , or possibly the treatment of nuclear cor- 
rections in deuterium, need to be reconsidered. 

The ratio of neutron to proton structure functions, 
F2/F2, is shown in Fig. [3] as a function of x* for various 
W* cuts (W* > 1.4, 1.6 and 1.8 GeV), and compared 
with the ratio from the recent CJ global PDF fit Q at 
matching kinematics. The range for the global fit arises 
from experimental and PDF fit uncertainties, as well as 
from uncertainties in the treatment of nuclear corrections 
in the analysis of inclusive data, which increase dra- 
matically at high x 0, Where the kinematics overlap, 
the data for the W* > 1.8 GeV cut are in good agreement 
with the global PDF fit for 0.3 < x* < 0.6 (the data at 
the lowest x* values arc outside of the range of validity 
of the global fit, which is restricted to Q 2 > 1.69 GeV 2 ). 
Note that a bump in F"£ / F^ appears near x* — 0.65 when 
relaxing the W* cut from 1.8 to 1.6 or 1.4 GeV, which 
likely indicates that a resonance in this region is signif- 
icantly enhanced in the neutron relative to the inelastic 
F$/F$ background. 

In summary, we have presented results on the first mea- 
surement of the neutron F£ structure function using the 



FIG. 3: (color online). Ratio F-f/Ff versus x* for various 
lower limits on W* . All data are from the 5.262 GeV beam 
energy. The error bars are statistical, with the total (cor- 
related and uncorrelated) systematic uncertainties indicated 
by the band along the abscissa. This band does not include 
the overall 3% normalization uncertainty or the 3% spectator 
approximation uncertainty. The data are compared with the 
recent parametrization from the CJ global analysis [f|, with 
the upper and lower uncertainty limits indicated by the solid 
lines. The inset shows the average Q 2 as a function of x* for 
each W cut. For these data a a is in the range 1.0-1.2. The 
arrow indicates the point at which the data are normalized 
to the CJ value. A single normalization constant Ivip was 
used for all data. The resonance region (W* < 2 GeV) cor- 
responds to x* > 0.4, 0.5 and 0.6 for square, diamond, and 
circle points, respectively. 



spectator tagging technique, where the selection of low- 
momentum protons at backward angles ensures scatter- 
ing from a nearly on-shell neutron in the deuteron. We 
identify well-defined neutron resonance spectra in each 
of the three prominent nucleon-resonance regions, which 
broadly agree with earlier modcl-dcpcndcnt extractions 
from inclusive deuteron and proton data but systemati- 
cally disagree in the details. The new, high-precision data 
will be useful in constraining models and paramctriza- 
tions of the neutron structure in the resonance region 
and beyond, and allow direct tests of quark-hadron du- 
ality in the neutron 0, [l(| ■ These will be the subjects of 
future publications. 

When combined with previous F^jF^ measurements, 
the new F^jF^ BoNuS data are used to reconstruct 
the ratio of neutron to proton F^ ' / F% structure func- 
tions up to x* » 0.6 in DIS kinematics, and up to 
x* w 0.8 in the resonance region, with little uncertainty 
due to nuclear effects. The results for the more strin- 
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gent W* > 1.8 GeV cuts agree well with the shape of 
recent global PDF fits [1, Q in regions where the kine- 
matics overlap, 0.3 < x* < 0.6, but show clear resonant 
structure at large x* for lower-IT* cuts. The precision of 
the new data, particularly in the DIS region, will be im- 
portant in reducing uncertainties in global PDF analyses 
13, 13 1 and extensions of the BoNuS experiment with the 
future 12 GeV Jefferson Lab will provide even stronger 

constraints on PDFs up to x « 0.8 [27 1. 
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